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Abstract 



The contribution of D waves to physical observables for neutral pion photopro- 
duction in the near threshold region is studied. Heavy Baryon Chiral Perturbation 
Theory to one loop, and up to 0(q 4 ), is used to account for the S and P waves, 
while D waves are added in an almost model-independent way using standard Born 
terms and vector mesons. It is found that the inclusion of D waves is necessary to 
extract the Eq + multipole reliably from present and forthcoming data and to assess 
the low-energy constants of Chiral Perturbation Theory. Arguments are presented 
C\| \ demonstrating that F-wave contributions are negligible in the near-threshold region. 



cn 

multipoles 

ON 
O 



Key words: pion photoproduction, chiral perturbation theory, electromagnetic 
multipoles 

PACS: 12.39.Fe, 13.60.Le, 25.20.Lj 



Due to the spontaneous breaking of chiral symmetry in Quantum Chromo 
Dynamics (QCD) the 7r meson appears as a pseudoscalar Nambu-Goldstone 
Boson pp. The dynamic consequences are that the S-wave amplitude for the 
7 A — > ir°N reaction is small in the threshold region, since it vanishes in the 
chiral limit, i.e. when the light quark masses are set equal to zero [2.3.4J. An 
additional consequence is that the P-wave amplitude is large and leads to the 
A resonance at intermediate energies [5]. Accordingly the photoproduction of 
neutral pions differs from the general pattern for hadronic reactions where the 
S wave dominates close to threshold and then, as the energy increases, the 
higher angular momentum states (P, D, . . . ) start to become important. By 
contrast, for the 7 A —* tt°N reaction the S- and P-wave contributions are 
comparable even very close to threshold [6]. 
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The purpose of this Letter is to show for the first time that the D waves also 
play a significant role in the near-threshold region (up to a photon energy 
~165 MeV) and to this end several observables have been extended beyond 
the traditional S+P wave limit to include the D-wave contributions. Fits have 
been performed for the latest and most accurate 7p — > ir p data in the near- 
threshold region [7J, in contrast to all previous analyses for which the data have 
been described using only S and P waves [2,3,4,7,6J. It will be shown that this 
affects the value of the extracted S-wave amplitude (E 0+ ), particularly its 
energy dependence. In contrast the P-wave multipoles are not affected by the 
inclusion of the D waves. 

Specifically, we explore the impact of D-wave contributions in the near-threshold 
region using as starting point Heavy Baryon Chiral Perturbation Theory 
(HBCHPT) [3.4J. We compute the S and P waves employing HBCHPT to 
one loop, and up to 0(q 4 ), and add the D waves using standard Born terms 
(equivalent to the Born contribution to HBCHPT) and vector mesons [8]. The 
Low Energy Constants (LECs) of HBCHPT are assessed through fits to the 
latest experimental data from Mainz [7J and the multipoles extracted. 

Let us begin by summarizing the basic structure of the observables and mul- 
tipoles. It is always possible to expand the differential cross section in terms 
of Legendre polynomials 

a T (9) = q - [T + T{P X (9) + T 2 V 2 (9) + . . . ] , (1) 

where and g 7 are the pion and photon momenta in the center-of-mass 
respectively and the Ts depend on the photon energy. Following the develop- 
ments in [9] we have worked out the full multipolar expansion of the Ts up 
to D waves. This expansion and a full analysis of near-threshold neutral pion 
photoproduction will be detailed in future work |10j . If we take into account 
only S and P waves, then only T , T\, and T 2 coefficients contribute. If we 
add D waves, two more quantities appear, T 3 and T 4 . So the first place to 
look for D waves consists in checking to see if there is room for the appear- 
ance of these new terms. The currently available experimental data [7j can be 
described quite well using only T , Ti, and T 2 , and no T 3 or T 4 contribution 
appears to be required at present. Hence, any contribution of D waves to o T 
should appear as a modification of To, T\ or T 2 . On the other hand, To and 
T 2 are dominated by diagonal terms involving the multipoles, namely |Mi + | 2 , 
|M!_| 2 , and |i?o+| 2 , and thus any interference with D waves would be negligi- 
ble compared to the leading terms. On the other hand, T\ is entirely due to 
multipole interferences: T\ up to D waves can be written: 

T l = 2Re[P?E 0+ ]+5T 1 , (2) 

where P\ = 3Ei + + M l+ — Ml_, and 5T\ stands for the D-wave/P-wave inter- 
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The P waves interfere with the D waves, enhancing their influence in the 
observable and possibly compromising any multipolar extraction from data. 
Note that there is no interference contribution to 5T\ between Eq + and the 
higher partial waves. 

To account for the S and P waves, the best available theoretical framework 
to study pion photoproduction in the near threshold region is HBCHPT. In 
[3lPf] the S and P multipoles to one loop and up to 0(q 4 ) are provided and we 
take these as our starting point. Six LECs appear, and five have been fitted 
to the data [7J : a\ and a 2 associated with the E 0+ counter-term, b p associated 
with the P 3 = 2M 1+ + M x _ multipole, together with £1 and £ 2 associated with 
P\ and P 2 = 3Pi + — M 1+ + M x _, respectively. The c 4 LEC associated with 
Pi, P 2 , and P3 has been taken from [11] where it was determined from pion- 
nucleon scattering inside the Mandelstam triangle. The corresponding details 
and fixed parameters, such as the pion-nucleon coupling g n ^, can be found 
in [3"P] . We obtain the D waves from the customary Born terms and vector 
meson exchange (u and p) [B] . For the uj and p parameters we have used those 
given by the dispersion analysis of the form factors in [12] . The vector-meson 
correction is very small and the inclusion of D waves in this way is almost 
equivalent to the zeroth order in HBCHPT. 

Within this framework we have performed two fits to the latest experimental 
data from Mainz [7] (171 differential cross sections and 7 photon asymmetries, 
spanning the energy range from threshold up to 166 MeV), one including 
only S and P waves (SP model) that was already done in [1], and another 
where we have added the D waves (SPD model). Both fits have the same 
number of parameters, which are the named LECs in HBCHPT. We obtain 
x 2 /dof = 1.23 for the SP model and X*/dof = 1.25 for the SPD model. 
The 70% (90%) confidence level is set to x 2 /dof = 1.27(1.29) and x 2 /dof = 
1.28(1.30), respectively. The P waves and their related LECs prove to be 
very stable against the inclusion of D waves, and indeed they are almost the 
same for both the SP and SPD fits and those given in [I], with differences 
lying typically below the 1% level. Of special interest is P\ which has £1 as 
its sole free LEC and is the one that could compromise the Eq + extraction. 
Both the SP and SPD models provide equally good fits to data from [7J. The 
P-wave parameters b p , £1, and £ 2 happen to be larger than those expected 
from resonance saturation. For an extensive analysis of this issue we refer the 
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Fig. 1. (Color online) Differential cross section at E~ = 154.2 MeV. The experi- 
mental points are from [TJ. The errors are statistical and do not include the 5% 
systematic error. We display results using the following models, SP (dotted), SPD 
(solid), and SP with D waves added without refitting the data (dashed). The SP 
and SPD curves mostly overlap. 
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Fig. 2. (Color online) T\ extracted from data [7J compared with results using the 
following models, SP (dotted), SPD (solid), and SP with D waves added without 
refitting the data (dashed). When theory is fitted to data including the D waves, 
the E$ + multipole is modified to compensate for the D-wave contribution. The SP 
and SPD curves mostly overlap. 

reader to [I]. The parameters of our SP fit differ from the ones in [4] and 
we obtain a better x 2 /dof. This is due to our better optimization procedure. 
We use a genetic algorithm combined with a standard gradient based routine, 
which provides us with a very powerful optimization tool. The details on the 
algorithm and its advantages can be found in [13] . In Fig. [1] we provide an 
example of how well the fits compare with experimental data from [7J. In 
addition we also provide the effects of D waves when added to the SP model 
without refitting. 

To obtain the experimental T\ displayed in Fig. [2] we have fitted each energy 
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Fig. 3. (Color online) Extracted Eq + multipole. The solid line corresponds to the 
SPD model and the dotted line to the SP model. 

data set from [7] to Eq. (CQ) using the least-squares method and assuming that 
only Tq, T 1; and T 2 contribute. In this way, we can extract T , T 1; and T 2 in an 
almost model-independent way (the only assumption being that higher-order 
Ts are negligible, as is the case given the current measurements). Together 
with the data, in Fig. [2] we plot T\ for the SP and SPD fits. They both fit 
the same data from [7] and so they provide a very similar T 1; although the 
multipolar contributions differ. As an exercise we also plot the SP model 
adding the D waves to explore how much they contribute to the observable. 
The modification of T\ is meaningful, and hence the SPD model becomes 
necessary. As a consequence of the stability of the P waves, the £"0+ has to be 
modified to compensate the D waves which pull T\ to more negative values. 
The curves show the influence of the unitary cusp in Re_Eo+, and are consistent 
with the data. To determine the magnitude of the cusp more accurately will 
require experiments with polarized targets |18j . 

When we analyze the fits, we find that P-wave multipoles are very stable 
against the inclusion of the D waves. However, E 0+ is significantly changed 
through a modification of the a± and a 2 parameters, i.e. via their impact on 
the counter-term (see below). Hence, the absence of D waves in the analysis 
compromises the extraction of the Eq + multipole and its associated LECs. In 
Fig. [3] we compare the two extractions of the E 0+ multipole. We can see that 
the slope and the unitary cusp are dramatically changed by the inclusion of 
the D waves, pointing out that D-wave contribution is meaningful even at the 
tt + threshold. 

Regarding which D-wave interferences play the major role in the modification 
of Ti, we have found that these are the interferences firstly of M 2+ and secondly 
of M 2 _ with the dominant multipole Mi+. At this point one might wonder 
if F waves could also make an important contribution to T\ in this energy 
region through their interference with P waves. The answer is no because, due 
to symmetry considerations, there is no interference of F and P waves in T]_. 
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Fig. 4. (Color online) Correlation plot for a\ and a 2 given the P-wave parameters. 
We display both the SP (dotted) and SPD (solid) models. The cross is the value of 
a\ and a 2 given by Bernard et al. [4j. For each model, the line is the 70% confidence 
level. Clearly seen are the strong correlation between the two parameters and the 
large uncertainties in their values. 



Indeed, F waves only interfere with D and G waves, guaranteeing a negligible 
contribution to T x in the near-threshold region. The argument is as follows: 
first, for T\ only partial waves with opposite parity can interfere (see Eq. ); 
second, these partial waves can only differ by one unit of angular momentum, 
unlike Tj with % > 1. Thus, one obtains the structure T\ = SxP + PxD + 
DxF + FxG + GxH + ... 

If we compare the SP and SPD fits for E 0+ we find that only two parameters 
can be modified, and those are ax and a 2 , the LECs that appear in the E 0+ 
counter-term: 

Eq + = ecuuMl + ea 2 u 3 , (4) 

where uo is the pion energy in the center-of-mass; these are significantly changed 
from one fit to another. In this analysis one has to keep in mind that the 
E 0+ multipole is a problematic one in HBCHPT due to its slow convergence 
[2113114115] . In principle these can be computed within Lattice QCD, and hence 
the extraction from data through HBCHPT becomes very important for the 
correct understanding of effective field theory. 

As noticed by Bernard, Kaiser, and MeiBner [3], a\ and a 2 are two highly 
correlated LECs, and their sum is a better quantity to use when analyzing 
data. In Fig. HJwe display the correlation plot for a\ and a 2 , where we can see 
that the two parameters are strongly correlated and that the error bars are 
extremely large. 

However, one can transform the counter-term from (ai,a 2 ) LECs to (a + = 
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In the near threshold region 
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Hence, the leading order for the counter-term is the one associated with a+, 
namely the parameter that can be well established by fitting data. At the 70% 
confidence level we obtain for the SP and SPD models, a + = 6. 481^93 GeV -4 
and a + = 6.57to;g2 GeV -4 , respectively 

In Fig. [5] we show the correlation plot for a+ and a_, where we can see how the 
a + can be well determined from the fits and is not affected by the inclusion 
of the D waves. On the other hand, a_ is strongly affected by the inclusion 
of D waves, changing completely its value and making the next order LEC 
a 2 = ( a + — a -)/2 very difficult to assess, as is shown in Figs. H]and0 

Another approach to extract the LECs is the one employed in p3], where the 
HBCHPT multipoles are fitted to the dispersion analysis of [15J to compute 
the Fubini-Furlan-Rossetti sum rule |16j . This approach has the shortcoming 
that it relies on the extraction made through a dispersion analysis with its 
own model dependencies and assumptions — it relies on the phenomenological 
model for pion photoproduction MAID03 [17J - - but has the advantages that 
HBCHPT has better convergence inside the Mandelstam triangle and that D 
waves are incorporated in the dispersion analysis. Actually, in [14J a truncated 
D-wave contribution was included in the fits and the LECs obtained are closer 
to the ones we find. We display the (a+, a_) pair from [H] in Fig. [5] as a plus 
sign. 

In summary, contrary to what is customarily claimed in the literature, the 
S+P approximation is not sufficient to obtain a complete description of the 
differential cross section in the near-threshold region and to extract the elec- 
tromagnetic multipoles reliably. The inclusion of D waves does not affect the 
extraction of the P-wave multipoles, but makes a big difference where the 
E 0+ extraction is concerned, especially in the energy region above the unitary 
cusp. The absence of D waves in the analysis compromises its extraction and 
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Fig. 5. (Color online) Correlation plot for a + and a_ given the P-wave parameters. 
We display both the SP (dotted) and SPD (solid) models. The cross is the value of 
a+ and a_ given by Bernard et al. |4j and the plus the value from Bernard et al. 
[14j . For each model, the inner line is the 70% confidence level and the outer one is 
the 90% line. Note that the different scales have been used for a+ and a_. 



the associated LECs of the Eq + counter-term. The current experimental in- 
formation is not yet accurate enough to pin down the a\ and a.2 LECs, but is 
capable of yielding a reasonable value for a + , allowing one to extract a reli- 
able Eq + at threshold to the order q A in HBCHPT. However, due to the slow 
convergence of the E 0+ multipole and the important effects of D waves and 
interferences one might wonder if the q A calculation is accurate enough and 
if a higher order computation may be necessary. The contribution of the D 
waves should be important in the interpretation of the new, more accurate, 
pion photoproduction data which was recently obtained by the Crystal Ball 
Collaboration at Mainz [19], as well as other future photo-pion production 
experiments with polarization degrees of freedom [T8f20"] . 
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